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1.  Introduction 

Lunar  physical  ephemerides,  tabulated  daily  including  libration,  can  be  found  in  Section  D  of  The  Astronom¬ 
ical  Almanac  (As A).  This  technical  note  describes  the  method  used  for  calculating  the  lunar  librations  and 
other  quantities  on  the  odd  pages  D7-D21  of  the  As  A  with  effect  from  the  2011  edition.  Prior  to  the  1985 
edition,  formulae  and  constants  for  the  physical  ephemeris  of  the  Moon  were  due  to  Hayn  (1907).  Beginning 
with  the  1985  edition,  librations  have  been  generated  using  the  analytical  theory  of  Eckhardt  (1981,  1982). 

Since  the  inclusion  of  a  rotational  ephemeris  for  the  Moon  given  in  the  form  of  three  Euler  angles  in  the 
JPL  lunar  ephemerides  starting  with  DE403/LE403,  Standish  et  al.  (1995),  it  has  been  a  long-standing 
goal  to  implement  these  ephemerides  in  the  calculation  of  lunar  librations  and  improve  the  quality  of  the 
tabulated  physical  ephemerides  of  the  Moon  in  Section  D.  Work  started  by  Hilton  (2004)  on  interpreting 
these  ephemerides  has  paved  the  way  to  the  work  presented  here  which  describes  the  calculation  of  librations 
and  related  quantities  as  they  have  been  implemented  starting  with  the  2011  edition  of  the  AsA.  This 
technical  note  has  been  written  with  the  practitioner  of  astronomical  calculations  in  mind  and  may  form  the 
basis  of  explanatory  material  on  lunar  librations  in  a  future  edition  of  The  Explanatory  Supplement  to  the 
Astronomical  Almanac  (ES). 

In  this  technical  note  the  basic  formulae  are  derived  for  the  lunar  librations  and  the  position  angle  of  the 
axis  of  rotation.  The  lunar  rotation  angles  in  the  JPL  ephemeris  must  be  transformed  by  a  series  of  rotations 
and  then  using  simple  vector  algebra  to  angular  quantities  which  can  be  substituted  into  the  basic  formulae. 
These  will  then  give  the  true  or  total  lunar  librations.  From  these  and  the  optical  (geometric)  librations,  the 
physical  librations  can  be  computed.  For  the  2011  edition,  the  lunar  rotation  angles  in  DE403/LE403  are 
used  due  to  the  availability  of  certain  ephemeris-specific  transformations  (Konopliv  et  al.  (2001)). 

An  algorithm  and  numerical  example,  using  the  same  computer  routines  used  for  the  AsA,  showing  all  the 
relevant  stages  in  the  computation  of  the  librations  and  related  quantities  is  given.  Appendix  C  contains  a 
short  description  of  the  calculation  of  the  position  angle  of  the  Moon’s  bright  limb  and  its  phase  or  fraction 
illuminated  -  geometrical  quantities  tabulated  on  the  odd  pages  D7-D21  of  Section  D. 

In  the  ES  (1961)  an  approximate  method  was  given  to  calculate  the  lunar  librations  but  with  no  details  of  its 
derivation.  In  Appendix  B  these  formulae  are  derived  based  on  the  method  in  Encke  (1843).  This  method 
was  used  primarily  before  the  advent  of  high-speed  electronic  computers  and  is  not  used  in  the  publications 
currently  to  calculate  the  lunar  librations.  It  is  included  here  for  the  interested  reader. 
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2.  Basic  formulae  for  the  lunar  libratioris  and- the  position  angle  of  the  axis 
2.1  The  lunar  librations 

The  mean  rotational  state  of  the  Moon  is  described  by  Cassini’s  empirical  laws,  which  state  that  the  descend¬ 
ing  node  of  the  Moon’s  equator  coincides  with  the  ascending  node  of  the  Moon’s  orbit  on  the  ecliptic;  the 
Moon’s  equator  mamtains  a  constant  inclination  to  the  ecliptic;  and  the  rotation  rate  is  such  that  on  average 
the  same  side  i$  always,  facing  the  .Earth.  Thus,  the, rotational  rate,  must  be  equal  to  the  rate,  of  motion  of 
the  Moon’s  mean  longitude.  The  actual  rotation  state  has  small  periodic  variations  from  this  mean  state 
caused  by  dynamical  perturbations,  and  these  cause  the  physical  librations  of  the  Moon’s  orientation.  In 
addition  there  are  the  much  larger  optical  librations  in,  its  orientation  as  seen  from  the  Earth,  which  are 
due  to  variations  in  the  rate  of  the  Moon’s  orbital,  motion,  and  to  the  inclination  of  the  Moon’s  equator  to 
its  orbital  plane  (see  Appendix  A  for  an  estimation  of  the  magnitude  of  these  librations).  For  a  complete 
definition  of  the  rotational  state  of  the  Moon,  a  prime  meridian  must  also  be  specified,  and  this  was  originally 
chosen  to  be  the  mean  central  meridian' of  "the  side  facing  the  Earth.  Its  direction  in  space  will  thus  differ 
by  180°  from  the  mean  longitude  of  the  Moon.  The  situation  is  illustrated  in  Figure  1. 


equator  on  the  ecliptic. 

The  ecliptic  longitude  and  latitude  of  the  Moon  axe  A, '  0  respectively  and  so  the  sub-Earth  point  M  has 
longitude  and  latitude  A  +  180°,  —  (3.  Lm  is  the  mean  longitude  of  the  Moon  and  Q  is  the  ascending  node. 
The  inclination  of  the  ecliptic  to  the  mean  lunar  equator  is  I.  The  librations  in  longitude  and  latitude  axe 
denoted  by  l  and  b  respectively. 

Formulae  for  computing"  the  optical  librations  can  be  derived  by  relating  to  each  other  two  expressions  for 
the  vector  from  the  centre  of  the  Moon  towards  the  sub-Earth  point  M;  one  of  them  referred  to  the  ecliptic 
frame  and  the  other  to  the' lunar  equatorial  frame.  Let  O  be  the  centre  of  the-selenOcentric  sphere.  In  Figure 
1  consider  firstly  point  M  referred  to  right-handed  axes  Oxyz  in  which  Ox  is  in  direction  OS,  Oy  in  the 
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ecliptic  and  Oz  in  direction  OC.  Coordinates  of  M  in  this  system  axe 


'cos(— 0)  cos(A  +  180°  —  12) 
cos(—0)  sin(A  +  180°  -  12) 
„  sin(— /?) 


(1) 


Consider  next,  point  M  referred  to  axes  Ox<y>i'  in  which  Ox>  is  in  direction  OS,  Oy>  in  the  mean-equator  of 
the  Moon  and  ;Oz'  in  the  direction  OL.  The  codrdinates  of  M  referred  to  these  axes  are 


'  cos  b  cos(Lm + 180°  —  12  + 1)  ’ 
cos  6sin(jLjtf  -h  186°  —  12-  + 1) 

, ' '  ■  sin-6  •  * 


(2) 


The  vectors  in  equations  (1)  and  (2)  are  related  through  the  rotation  matrix  . 

I1  0  0  V 

Ri(— I)— -I  0  -cos I  —sin/  1 
\  0  sin  I  cos'/  , 


(3) 


We  note  here  in  general  that  rotations  in  a  right-handed  set  of  axes  with  origin  O  about  the  Ox,Oy  and  Oz 
axes  through  an  arbitrary  angle  0  are  obtained  by  the  rotation  matrices  Ri,  R2  and  R3  respectively,  where 


0  O' 
cos  9  sin  9 
—  sin  9  cos  6  / 

cos  9  0  —  sin  9 ' 
0  10 
sind  0  cos  9  / 

cos  8  sin  8  0 ' 

—  sin  9  cos  9  0 

0  0  1, 


(4) 

(5) 

(6) 


Prom  equations. (1),  (2)  and  (3). we  have 

/  cos  b  cos(L m  +  180“  —  1?  +/)  \ 

I  cos  b  sin {LM  + 180°  +  0  +4) ..  I  = 
\  sin  6  •  J 

which  can  be  written  as 

(  —  cos  6  cos(/'M  +-i  — /2)\ 

—  cos  6sin(LM  +1  —  0)  J  = 
sin  6  J 


/  cos(—0)  cos(A.+T80°  —  12) 
i  cos(^3);sm(A  +  l'80°  -  1?) 
\  '  sin(-/3) 


1  0  0 
0  cos  I  —  sin  I 
.  0  sin  /  cos  / 


-cos/?  cos  (A  —  12) ' 
-  cos  0  sin(A  —  12) 
•—  sin  0 


(7) 


In  the  lunar  ephemeris  A  is  referred  to  the  true  equinox  of  date,  but  Lm  and  12  are  referred  to  the  mean 
equinox.  The  quantity  A  —  12  in  equation  (7)  must  be  replaced  by  A  —  (12  +  N),  where  N  is  the  nutation  in 
longitude.  We  can  now  write'  (7)  as  the  three  equations 

cos6cos(i  +  Lm  —  12)  =  cos/?cos(A  —  12  —  N)  •  (8) 

cos6sin(Z  +  Lm.  —  12)-=  cos/cos/3si.n(A 12  — -JV)  —  sin/sin/J  (9) 

sin  b  sin  /  cos  0  sin(A  —  12  —  Ar)  —  cos  I  sin  0.  ■  (10) 

Equations  (8),  (9)  and  (10)  axe  rigorous  formulae  for  the  computation  of  the  optical  librations  l  and  6,  from 
the  values  of  I,  12  and  Lm,  which  describe  the  mean  rotational  state;  of.the  Moon,;  and  from  A  and  0,  the 
ecliptic  coordinates  of  the  Moon.  However,  as  is  explained  in  Section  3;  .if  /,  12.  and  Lm  are  substituted 
by  modified  quantities  that  include  the  effects  of  the  dynamical  perturbations  of  the.  Moon’s  rotation;  then 
these  rigorous  formulae  will  give  the  values  of  the  total  librations. 
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An  approximate  method  to  compute  the  libfations  l,  b  was  given  in  the  ES  (1961)  p.319  based  on  formulae 
introduced  by  Encke  (1843).  It  is  derived  using  the  quantities  B  and  l*  shown  in  Figure  1  and  was  used 
before  the  advent  of  fast  electronic  computers.  A  statement  of  the  method  and  an  outline  of  its  derivation 
is  given  in  Appendix  B  for  the  interested  reader. 

2.2  The  position  angle  of  the  axis 

The  position  angle  of  theaxis  of  rotation  is  the  angle  that  the  lunar  meridian  through  the  apparent  central 
point  of  the  disk  towards  the  north  lunar  pole  forms  with  the  celestial  meridian  through  the  central  point, 
measured  eastwards  from  the  celestial  north  point  of  the  disk. 

In  determining  expressions  for  the  position  angle  we  use  the  elements  of  the  mean  lunar  equator  referred  to 
the  Earth  equator. 


Figure  2:  Elements  for  Moon’s  equator.  , 

These  are  defined  as:\  :  -  . 

"  r*.  l  • 

,  '  i  =  the  inclination  of  the  mean'  equator  of  the  Moon 
/  to  the  true  equator. of  the  Earth;'  . 

s  A  =  the  arc  of  the  mean' equator  of  the  Moon  from  its 

ascending,  node  bn  the  true  equator  of  the  Earth 
to;  it's' ascending,  node  on  the  ecliptic  of  date; 

•  Q'  =  the  arc. of  the  true  equator  of  the  Earth  from -the 

true  equinox  of  date  to  the  ascending  node  of  the 
mean  equator  of  the  Moon  on  the  true  equator 

j.’v 

’■  '  of  the  Earth. 


The  ascending  node  of  the  mean  lunar  equator  on. the  ecliptic  is  at  the  descending  node  of  the  mean  lunar 
orbit  so  TS'  —  C2±  1809.  e  is  the  true  obliquity  and  the  node  is  referred  to  the  true  equinox  by  increasing 
Q  by  the  .nutation  in  longitude  N.  From  the  spherical  triangle  TS’tj'  in  Figure  2  the  elements  can  be  found 
from  the  formulae  ■  ■•■  ■  ■■ 


■  sin^sini  —  —  sine-sin (12  +  N)  ■  (11) 

■  'bos  2lsini  =  sin/cose  —  cos  Isin£cos(I2  4-  N)  •  .  (12) 

cost  =  cos  loose  ■  t~sin  I  sin£cos(£?--F  Ar)  (13) 

sin  Q'  sinf  —  —  sin  I  sin(I2  +  N)  (14) 

cos  12'  sini  =  cos 7 sins  —  sin/cosecos(I2  +  N).  '  (15) 


In  Figure  3  the  position  angle  C'  of  the' axis  is, shown  on  . the  selenocentric  sphere.  The  geocentric  right 
ascension  and  declination  of  the  Moon  are  a,  5  and  so  the  right  ascension  and  declination  of  the  sub-Earth 
point  are  a  +  180°,  ~<5.  The  descending  node  of  the  lunar  equator  on  the  ecliptic  is  denoted  by  point  S. 
From  the  definition  of  A  and  Q'  the  arcs  SQ  =  360°  —  A  and  TQ  =  Q'  -I- 180° .  In  the  spherical  triangle 


4 


NLM  all  angles  and  sides  axe  known  except  angles  LNM  and  MLN.  Noting  that  XQ  —  90°  and  YQ  =  90° 
these  are  found  as  follows: 

LNM  =  180°  -YP 

=  180 °-{YP-TY) 

—  180°  -  TP  +  (YQ  -  YQ)  -  ' 

=  180°  ~{ot  +  180°)  +  (12'  +  180°) -90° 

■■  =  n'  -  a  +  90°  .......  (16) 

" .  MLN  =  XT  " 

—  XQ  —  TQ  *'■ 

'  ■  ’  =XQ~(SQ-  ST) 

=  90°  -  (360°  -  A)  +  {Lm  - 12  + 180°  4- 1) 


:=A  +  LM  +  l-  n- 90°.  (17) 

The  position  angle  (C)  can  be  found  from  either  of  the  two  sets  of  formulae 

cos  b  sin  C'  =  — '•sini  cos  (12'  —  a)  (18) 

cosh  cos  C'  =  cos  5  cos  i  —  sin  5  sinz  sin(f2'  —  a)  (19) 

or 

cos  <5 sin  C"  =  sinicos(LM-  —  Q  +  A  ±1)  (20) 

■cos fipps C'.~cgsicgsb  —  sin  z  sin  6  sin  (£71* r—  fi  +  A  4- l).  (21) 


Figure  3:  The  selenocentric  sphere:  showing  the  relationships  between  the  sub-Earth  point  M,  the 
mean  lunar  equator  and  the  Earth’s  equator. 
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3.  Implementation  of  JPL  lunar  rotation  angles 

The  JPL  DE403/LE403  ephemeris,  commonly  known  as  DE403,  includes  an  ephemeris  for  the  rotation  of 
the  Moon.  The  problem  is  to  transform  these  rotation  angles  into  variables  equivalent  to  the  quantities  Q , 
I,  Lm  and  hence  use  equations  (8),  (9),  (10),  and  (18),  (19)  (or  (20),  (21))  with  these  newly  determined 
quantities  to  compute  librations  It,  fey  and  position  angle  C't-  We  refer  to  It,  as  the  true  or  total  lunar 
librations.  ‘ 

The  Euler  angles  describing  the  rotation  . of  the  Moon,  <p,  9  and  ip  are  determined  in  Newhall  and  Williams 
(1997).  These  angles  are  defined  relative  to  the  .  ICRS  Earth  equator  and  equinox.  They  describe  the 
orientation  of  the  principal  axes  of  inertia  of  the  Moon  (the  PA  system),  sometimes  called  the  axes  of  figure 
system.  We  show  how  they  can  be  transformed  to  give  new  Euler  angles  <pc,  and  ipc,  which  are  defined 
relative  to  the  ecliptic  and  equinox  reference  frame  .of  date.  These  transformed  angles  are  used  to  describe 
the  orientation  of  a  slightly  different  lunar  axis  system,  which  has  one  axis  towards  the  mean  Earth  direction, 
'  and  another 'along  th'e  rotation  aid's -(the  ME'Sy'steifi), 'sometimes  called  the  mean  Earth/ rotation  axis  system. 
From  the'se  'quantities -we  can  cttmpute  new  S2, T  and  Lm--  ■  ■  ■  ' 

The  Euler  angles  6,  r/>,and.  new  .set  pci  6c,.ipc  are,  shown  in  Figures  4a,  4b  respectively.  .These  are  given 
in  Newhall  and.  Williams  (1997)  and gre, included  in  this  .technical .note  for  convenience.. 


Figure  4a.:  Equatorial  ref erence  frame  showing  the  Euler  angles  -<p,  9,  ip,  used  to  describe  the  lunar 
principal  axis  (PA )  system.  The  value  of  <p  shown  in  this  diagram  is  negative. 


Figure  4b:  Ecliptic  reference  frame  showing  the_  transformed  Euler  angles  <pc,  6c,  'Pc,  nsed  to 
describe  the  orientation  of  the  lunar  ME  system. 

The  angles  are  defined  as: 

4>  =  the  angle  along  the  ICRS  equator,  from  the  ICRS  X-axis  to  the  ascending 
node  of  the  lunar  equator; 

9  —  the  inclination  of  the  lunar  equator  to  the  ICRS  equator; 
ip  =  the  angle  along  the  lunar  equator  from  the  node  to  the  lunar  prime  meridian; 
pc  ~  the  angle  from  the  equinox  of  date  to  the  descending  node  of  the  lunar 
equator  on  the  ecliptic  of  date; 
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6c  =  the  inclination  of  the  lunar  equator  to  the  ecliptic  of  date; 
ipc  =  the  angle  along  the  lunar  equator  from  its  descending  node  on  the  ecliptic 
to  the  lunar  prime  meridian. 


From  Seidelmann  et  al.  (2007)  a  vector  p  in  the' PA  system  can  be!  transformed' to  a  vector  q  in  the  ME 
system  by  applying  three  small  rotations.  This  expression  (Konopliv  et  al.  (2001))  is 

q  =  Ri.(— 0^T4&2)  R2(— 79<'0768)  R3  (-63^8986)  p  ■  (22) 

where  Ri ,  R2 ,  R3  are  given  in  equations  (4) y  (5)  ,-  (6)  ■  resp  ectively.  It'  must  be  noted  that  the  numerical 
values  for  the  rotations  in  equation  (22)  are  specific' to  DE403  and’  are  different  for  other  ephemerides. 

The  JPL  lunar  libration.ephem.eris  gives  the  orientation  of.  the  PA  .syst.em  relative  to  the  ICRS.  So  if  we 
know  the  components  of  any  vector  in  the  PA. system,  then  the.  JPL  libration  angles  enable  ns  to  refer  it  to 
the  ICRS.  However  the  commonly- used  reference  system  for  lunar  cartography  is  the  ME  system.  Using  the 
inverse  of  equation  (22)  we  can  convert  any  vector  in  the  ME  system  into' the  PA  system;  and  then  we  can 
use  the  JPL  libration  angles  to  convert  to  the  ICRS:  Finally  we  apply  frame' bias,  precession  and  nutation, 
and  convert  to  the  ecliptic  reference  frame  of  date.  We  apply  this  process  to  vectors  in  the  ME  system  that 
define  the  lunar' pole  and  tire  lunar  prime  meridian,  and  use  the  converted  vectors  to  calculate  lunar  libration 
angles  for  the  ME  systeih  in  the  ecliptic-of-date  frame.  ..  .  " 

Let  rsel  be  a  vector  in  the  ME  system- The  vector  ri  in  the  PA  system  is  found  from  the  inverse  of 
equation  (22)  . * -  '  ■ '■ . 

ri  =  R3_1  (-63^8986)  Ra-1  (-79!,07fi8)'R1_1  (-0^1462)  rseI 

=  R3(63"8986)  Ra  (79^0768)  Rj  (0!'1462)  rseI.  (23) 

In  the  ICRS  equator- system  (see.  Figure '4(a))  this  vector  becomes 

r2.='R.3(-0)  Ri(-0)  R3(-V0  (24) 

Finally,  in  the  ecliptic  of  date  system  this  vector  becomes 

rdate=R,(e)NPBr2  ...  .  (25) 


where  B,  P  and  N  are  the  frame  bias,  precession  and  nutation  matrices  respectively,  and  Rj(e)  is  the 
rotation  to  the  true  ecliptic,  of.  date  reference  frame.  Note  that  B,  the  frame  bias  matrix,  is  included  for 
completeness,  as  its  effect  is  well  under  O-'l. 

Now  apply  rotations  in  equations  (23),  (24),  (25)  to  r5®1  =  (1,0,0)  and  then  to  r5el  =  (0,0, 1)  and  call  the 
resulting  vectors  x-date  and  zdate.respectivcly.ythey  are.  with'  respect  to  the' ecliptic  of  date  system. 

In  the  ecliptic  of  date  system  let  k  b'e  unit  vectors  along  the  Ox,  Oy  and  Oz  axes  respectively.  Define 
the  unit  vector  f2  to  be  in -the  direction  0  to  the  descending  node  of  the  lunar  equator  on  the  ecliptic  (S). 
We  have 


zdatexk 
|  zdatexk  | 


(26) 


The  angles  <j>c,  @c  and  ipc  are  found,  using  Figure  4(b),  from  the  following  formulae 

cosific  =  i.fi  (27) 

sind’c  =  j-O  .-  ,  ,  (2S) 

cos  6c  =  k.zdate  (29) 

cos  tpc  =  f2.xdate  (30) 

'  silled  =  (zdatcxf2).xdate.  (31) 


From  Newhall  and  Williams  (1997) 

4>c  12  -1-  c 

0c  =  I  +  p  ■ .  - .  .  •  ■ 

ipc  =  t  —  a  +  F  4- 180°.-  v. 
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(32) 

(33) 

(34) 


Since  (see  ES  (1961)  p.107) 


F  —  Lm  —  12 


equation  (34)  becomes  using  (35) 


ipc  =  t  ~  cr  +  Lm  —  12  +  180°. 


Substituting  now  for  Q  from  equation  (32)  into  (36)  we  have 

tpc  =  r  ~  <t  +  Lm  —  4>c  +  o’  + 180° 
'  tfc  +  4>c  —  180°  =  Lm  4-  t. 


(35) 

(36) 


(37) 


Prom* equations  (32),  (33)  and  (37)  the  changes  in  12, 1  and  Lm  viz-,  oypand  t  respectively  are  as  a  result 
of  the  physical  libration.  To  obtain  the  total  or  true  librations  lT,  br  we  must  therefore  substitute-^.  @C 
and  4>c  +  <Pc  —  lSO.Vfor  12, 1  and  Lm  respectively  in -equations  (8),-(9)  and  (.10).  We  note- since  the  values 
for  12  and  Lm  that  aresubstituted-into  equations  (8),  -(9)  and  (10)  are  referred  to  the  true:  equinox  we  must 
set  N  =  0.  The  position  angle  C't  for  the  total  librations  is  then  found  using  equations  (18),  (19)  or  (20), 
(21).  •  •  •  -  1  " 

We  then  can  calculate  the  physical  librations  dip,  Sbp  arid  5C'P  from 


dip  =  It  —  lb  (38) 

•  •;  ..  .  ,  ■■  SbP  —  bp  —  bo  .  . .  ,  ■  (39) 

SC'p  =  C'T  -  C'o  (40) 

where  the  optical  librations  lo,  bo  and  the. position  angle]  Q'o  are  as  computed  from  equations  in  section  2. 

The  heliocentric. ecliptic, longitude,  and.  latitude,  . fig,  of  the.Mpon  are  determined  by  calculating  the 

vectors  for  the  geocentric  ecliptic  positions  of  the  Sun  and  Moon  and  forming  the  heliocentric  vector  to  the 
Moon'll Y, Z) <sm -'artd'itS'length,  d$-Mi  "To  :-.V.  ■  .. 


'"  /cos'X'C6sf3\  ■  N  • 
d  1  sin  A  cos /3  ’  -  •  '■ 

si v,P  .  / 


(X  V  /  cos  As  cos  ps  V 

Y  =ds  sin  As  cos  0s  ] 
ZJS  .\~:<sm0s  ) 


dsM  — 


+:  %SM 


tan  A/f  - 


Ysm  ■  '■ 
XsM 


sin = 


Zsm 

dsM 


(41) 

(42) 

(43) 

(44) 

(45) 

(46) 


By  substituting  A#  for  A  and  ftp  for  /?  in  equations  (8)  and  (9),  the  'stenographic  longitude  of  the  Sun, 
Is,  can  be  determined  in  a  similar  manner  to  the  libration  in  longitude.  Similarly,  substitution  of  these 
quantities  into  equation  (10)  allows  the  determination  of  the  selenographic  latitude,  bs,  in  a  similar  manner 
to  the  libration  in  latitude.  The  selenographic  colongitude  of  the' Sun  is  simply  90°  —  is,  adjusted  to  lie  in 
the  range  0°  to  360°.  ,  .. 

4.  Numerical  Example 

The  purpose  of  this  numerical  example  is  to  calculate  the  quantities  on  the  odd  pages  D7-D21  of  the  AsA, 
most  of  which  involve  the  lunar  librations.  The  quantities  tabulated  are: 
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•  The  Earth’s  selenographic  longitude  (It)  and  latitude  (bp) ,  which  correspond  to  the  total  librations 
in  longitude  and  latitude,  respectively. 

•  The  position  angle  of  the  axis  of  rotation  (Cj.). 

•  The  physical  librations,  5lp,  6bP  and  the  difference  8C'P,  tabulated  in  thousandths  of  a  degree. 

•  The  Sun’s  selenographic  colongitude  (90°  —  l $)  and  latitude  (65). 

•  The  geometrical  quantities  fraction  illuminated  (/j)  and  the  position  angle  of  the  bright  limb  (PA p). 
4.1  Algorithm  for  calculating  the  librations  and  position  angle  of  the  axis 


The  method  used  involves  two  passes  through  a  common  process.  The  first  pass  calculates  the  optical 
librations,  which  are  geometrical  innature-and  adhere  to 'Cassini’s  laws.  The  second  pass- calculates  the 
total  librations  which  includes  the  adopted  rotational  ephemeris  of  the  Moon.  In  this  example,  the  JPL 
DE403  Ephemeris  is  used  (see  Section  3).  I  AU  Standards- of  Fundamental  Astronomy  routines  {SOFA)  are 
used  for  obtaining- the  frame  bias,  precession  and  nutation.  The  steps  are  given  below.  • 

Step  1  Obtain  the  apparent  positions  a,  5,  d,  and "a3,  5S,  ds  for  the  Moon  and  Sun,  respectively,  at  time  t. 

Step  2  Obtain  the  nutation  in  longitude  (Aqb)  and  obliquity. (Ae)  and  the  true  obliquity  (s)  for  time  t. 

Step  3  Determine  the  apparent  ecliptic  positions  A,  0,  and  A.,,  0Sl  by  rotating  the  apparent  equatorial 
coordinates  around  the  A-axis  by  the  angle  s. . 

Step  4  Calculate  the  light  time  correction  fdr  the' Moon,  r  =  d/c,  where  c  =  173-1446  3268  467  is  the  speed 
of  light  in  au/day.  ■  ' 

Step  5  At  time  t  —  t  form  the  true  inclination  Im  +.  As,  where  Ae  is  the  nutation  in  obliquity  and  Im  is 
the  Newhall  and  Williams  (1997)'  value’of  the  inclination.  ' 

Step  6  Also  at  time  t  .—'T  evaluate  the  fundamental  arguments  12  and  L'm  Using  Simon  et  al'(1994). 

Step  7  Using  equations  (8),  (9)  and  (10),  calculate  the  optical  (geometric)  librations  1  and.b  from.  \—C2  —  N, 
0,  I  and  Lm ,  where  N  =  Aip  is  the  nutation  in  longitude. 

Step  8  Next,  in  order  to  determine  C1,  the  position  angle  of  the  axis  of  rotation,  determine  Q'  from  equations 
(14)  and  (15),  and  i  from  equation  (13)-,  Thus  C  may  be  calculated  from  equations  (18)  and  (19). 

Step  9  Now  use  the  ephemeris  and  equations  (27)-(31)  to  determine  the  Euler  angles  <j>c,  9c  and  ipc-  This 
process  includes  the  transformations  from  the  system  pf  the  ephemeris  (JPL  DE403)  to  the  true  of 
date  system  equations  (23)-(25).  '  . ' 

Step  10  Then  substitute  4>c  for  Q,  6c  for  1  and  bf>c  +  <j>c  —  180°  for  Lm  and  repeat  steps  (7),  (8)  and  (9) 
resulting  this  time  in  the  total  librations  Ip,  bp  and  C'T. 

Step  11  Lastly  the  physical  librations,  the  differences  between  the  total  and  the  optical  librations  are  given 
by  equations  (38)- (40).  ■■ 

'  ' 

4.2  Numerical  Example  .._  _  . .. 

Calculate  the  quantities  described  in  this.  Technical -Note  for;  2011  June  01  at  OhTT,  approximately  21  hours 
before  the  instant  of  new  Moon. 

t  =  2455713-5000000  2011  June  01  0h  TT 


Step  1 ;  The  apparent  positions  of  the  Moon  and  Sun  have  been  calculated  using  methods  described  in 


Section  B  of  The  Astronomical  Almanac. 


Moon 

a  — 

57^364896851 

5  = 

22?200527037 

d  = 

0-0026441632  au 

Sun 

as.  = 

68°b64159796 

5s  =  ., 

211975380381 

ds  = 

1-0139593548  au 

Step  2:  Nutations  and  the  mean  and  true  obliquity  have  been  calculated  using  SOFA  routines. 


IV  =  Aip  = 
As  = 

£m  = 

£  =  £m  -H  Ae  = 


0^004500032 

—0^000366339 

23M37794624 

23°437428285 
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